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Abstract 

The production of D* + , D°, D + , and A+ charm hadrons and their antipar- 
ticles in ep scattering at HERA was measured with the ZEUS detector using 
an integrated luminosity of 79pb _1 . The measurement has been performed in 
the photoproduction regime with the exchanged-photon virtuality Q 2 < 1 GeV 2 
and for photon-proton centre-of-mass energies in the range 130 < W < 300 GeV. 
The charm hadrons were reconstructed in the range of transverse momentum 
Pt(D,A c ) > 3.8 GeV and pseudorapidity \rj(D,A c )\ < 1.6. The production cross 
sections were used to determine the ratio of neutral and charged D-meson pro- 
duction rates, R u /d, the strangeness-suppression factor, j s , and the fraction of 
charged D mesons produced in a vector state, P*. The measured R u /d and 7 S 
values agree with those obtained in deep inelastic scattering and in e + e~ an- 
nihilations. The measured P^ value is smaller than, but consistent with, the 
previous measurements. The fractions of c quarks hadronising as a particular 
charm hadron, f(c — * D,A C ), were derived in the given kinematic range. The 
measured open-charm fragmentation fractions are consistent with previous re- 
sults, although the measured f(c — > D* + ) is smaller and f(c — + A+) is larger 
than those obtained in e + e~ annihilations. These results generally support the 
hypothesis that fragmentation proceeds independently of the hard sub-process. 
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1 Introduction 



Charm quark production has been extensively studied at HERA using D* ± and Df 
mesons [1-5]. The data have been compared with theoretical predictions by assuming 
the universality of charm fragmentation and using the charm fragmentation characteris- 
tics obtained in e + e~ annihilation for the calculations of charm production in ep scat- 
tering. However, the charm production mechanisms are not the same in different colli- 
sions. In particular, cc pairs in e + e~ annihilation are produced dominantly in a colour- 
singlet state, which is not the case for ep scattering. Thus, it is important to test the 
charm-fragmentation universality by measuring the charm fragmentation characteristics 
at HERA. 

In this paper, the measurement of the production of the weakly decaying charm ground 
states, the D°, D + , pseudo-scalar mesons and the A+ baryon, is presented. The 
production of the charm vector meson D* + has also been studied. The antiparticles 
of these charm hadrons have been measured as well 1 . The measurement has been per- 
formed in ep scattering at HERA in the photoproduction regime with exchanged-photon 
virtuality, Q 2 , close to zero and for photon-proton centre-of-mass energies in the range 
130 < W < 300 GeV. The measured production cross sections have been used to deter- 
mine the ratio of neutral and charged D meson production rates, R u /d, the strangeness- 
suppression factor, 7 S , and the fraction of charged D mesons produced in a vector state, 
Py. The fractions of c quarks hadronising as a particular charm hadron, f(c — > D,A C ), 
have been calculated in the accepted kinematic range. The open-charm fragmentation 
fractions in photoproduction are reported here for the first time. The results have been 
compared with the previous HERA measurements of the charm fragmentation character- 
istics in photoproduction [4] and in deep inelastic scattering (DIS) with Q 2 > 2 GeV 2 [6]. 
To compare the results with those obtained in charm production in e + e~ annihilations, the 
/(c — > D, A c ) fractions compiled previously [7] have been updated using recent values [8] 
of the relevant branching ratios. 



2 Experimental set-up 

The analysis was performed with data taken by the ZEUS Collaboration from 1998 to 
2000. In this period, HERA collided electrons or positrons 2 with energy E e = 27.5 GeV 
and protons with energy E p = 920 GeV. The results are based on a sum of the e~p and 
e + p samples corresponding to a total integrated luminosity of 78.6 ± 1.7 pb~ . Due to 

1 Hereafter, charge conjugation is implied. 

2 From now on, the word "electron" is used as a generic term for electrons and positrons. 
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trigger considerations, D + and A+ production was measured using only the e + p sample 
corresponding to an integrated luminosity of 65.1 ± 1.5 pb~ . 

A detailed description of the ZEUS detector can be found elsewhere [9]. A brief outline 
of the components most relevant to this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [10], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organized in nine superlayers covering the 
polar-angle 3 region 15° < 6 < 164°. The transverse-momentum resolution for full-length 
tracks is a(p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. To estimate the 
energy loss per unit length, dE/dx, of particles in the CTD [4,11], the truncated mean of 
the anode-wire pulse heights was calculated, which removes the lowest 10% and at least the 
highest 30% depending on the number of saturated hits. The measured dE/dx values were 
normalised to the dE/dx peak position for tracks with momenta 0.3 < p < 0.4 GeV, the 
region of minimum ionisation for pions. Henceforth dE/dx is quoted in units of minimum 
ionising particles (mips). The resolution of the dE/dx measurement for full-length tracks 
is about 9%. 

The high-resolution uranium-scintillator calorimeter (CAL) [12] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/ \J~E for electrons 
and a(E)/E = 0.35/ \[E for hadrons, with E in GeV. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > e'yp, 
where the photon was measured with a lead-scintillator calorimeter [13] located at Z — 
-107 m. 

3 Event simulation 

Monte Carlo (MC) samples of charm and beauty events were produced with the Pythia 
6.156 [14], Rapgap 2.0818 [15] and Herwig 6.301 [16] event generators. The generation, 
based on leading-order matrix elements, includes direct photon processes, in which the 
photon couples as a point-like object in the hard scatter, and resolved photon processes, 

3 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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where the photon acts as a source of partons, one of which participates in the hard 
scattering process. Initial- and final-state parton showering is added to simulate higher- 
order processes. The CTEQ5L [17] and GRV LO [18] par ametrisat ions were used for 
the proton and photon structure functions, respectively. The charm and bottom quark 
masses were set to 1.5 GeV and 4.75 GeV, respectively. Events for all processes were 
generated in proportion to the predicted MC cross sections. The Lund string model [19] 
as implemented in Jetset [14] was used for hadronisation in Pythia and Rapgap. The 
Bowler modification [20] of the Lund symmetric fragmentation function [21] was used for 
the charm and bottom quark fragmentation. In Herwig, the cluster model [22] was used 
for hadronisation. The fraction of charged D mesons produced in a vector state was set 
to 0.6 for all MC samples. 

The Pythia and Rapgap generators were tuned to describe the photoproduction and DIS 
regimes, respectively. Consequently, the Pythia events, generated with Q 2 < 0.6 GeV 2 , 
were combined with the Rapgap events, generated with Q 2 > 0.6 GeV 2 . Diffractive 
events, characterised by a large rapidity gap between the proton at high rapidities and 
the centrally-produced hadronic system, were generated using the Rapgap generator in 
the diffractive mode and combined with the non- diffractive MC sample. The contribution 
of diffractive events was estimated by fitting the 77 max distribution 4 of the data with a linear 
combination of the non-diffractive and diffractive MC samples. The combined sample was 
used to evaluate the nominal acceptances. The Herwig MC sample, generated over the 
full range of Q 2 values, was used to estimate the model dependence of the acceptance 
corrections. 

To ensure a good description of the data, the transverse momenta, pt{D,A c ), and pseu- 
dorapidity, rj(D, A c ), distributions were reweighted for both combined Pythia+Rapgap 
and Herwig MC samples. The reweighting factors were tuned using a large D*^ sam- 
ple [23] . The effect of the reweighting on the measured fragmentation ratios and fractions 
was small; the reweighting uncertainty was included when estimating the model depen- 
dence of the acceptance corrections. 

The generated events were passed through a full simulation of the detector using Geant 
3.13 [24] and processed with the same reconstruction program as used for the data. 

4 Event selection 

A three-level trigger system was used to select events online [9,25]. The first- and second- 
level trigger used CAL and CTD data to select ep collisions and to reject beam-gas 

4 The quantity ?7max is denned as the pseudorapidity of the CAL energy deposit with the lowest polar 
angle and an energy above 400 MeV. 
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events. At the third level, where the full event information was available, at least one 
reconstructed charm-hadron candidate was required. The efficiency of the online charm- 
hadron reconstruction, determined relative to the efficiency of the offline reconstruction, 
was above 95%. 

Photoproduction events were selected by requiring that no scattered electron was identi- 
fied in the CAL [26]. The Jacquet-Blondel [27] estimator of W, W ]B = ^1E V (E — p z ), 
was used, where E — pz = — pz)i and the sum % runs over all final state energy-flow 
objects [28] produced from charged tracks, as measured in the CTD, and energy clusters 
measured in the CAL. After correcting for detector effects, the most important of which 
were energy losses in inactive material in front of the CAL and particle losses in the beam 
pipe [26,29], events were selected in the interval f30 < W < 300 GeV. The lower limit 
was set by the trigger requirements, while the upper limit was imposed to suppress re- 
maining DIS events with an unidentified scattered electron in the CAL [26]. Under these 
conditions, the photon virtuality lies below f GeV 2 . The median Q 2 value was estimated 
from a Monte Carlo simulation to be about 3 x f -4 GeV 2 . 

5 Reconstruction of charm hadrons 

The production of D* + , D°, D + , and A+ charm hadrons was measured in the range of 
transverse momentum pt(D,A c ) > 3.8 GeV and pseudorapidity \rj(D,A c )\ < 1.6. Charm 
hadrons were reconstructed using tracks measured in the CTD and assigned to the re- 
constructed event vertex. To ensure good momentum resolution, each track was re- 
quired to reach at least the third superlayer of the CTD. The combinatorial background 
was significantly reduced by requiring p T {D) / E^ >10 ° > 0.2 and pt{A c )/E^ >10 ° > 0.25 
for charm mesons and baryons, respectively. The transverse energy was calculated as 
£>0>io° _ j]j i6 i. >10 o(£'jsin^j), where the sum runs over all energy deposits in the CAL with 
the polar angle 9 above 10°. Further background reduction was achieved by imposing 
cuts on the transverse momenta and decay angles of the charm-hadron decay products. 
The cut values were tuned using MC simulation to enhance signal over background ratios 
while keeping acceptances high. 

The details of the reconstruction of the five charm-hadron samples are given in the next 
sub-sections. 

5.1 Reconstruction of D° mesons 

The D° mesons were reconstructed from the decay D° — > K~ir + . In each event, tracks 
with opposite charges and pt > 0.8 GeV were combined in pairs to form D° candidates. 
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The nominal kaon and pion masses were assumed in turn for each track and the pair 
invariant mass, M(Ktt), was calculated. The distribution of the cosine of the D° decay 
angle (defined as the angle 8*(K) between the kaon in the Kir rest frame and the Ktt 
line of flight in the laboratory frame) is flat, whereas the combinatorial background peaks 
in the forward and backward directions. To suppress the background, | cos9*(K)\ < 0.85 
was required. 

For selected D° candidates, a search was performed for a track that could be a "soft" pion 
(jTs) in a D* + — > D°irf decay. The soft pion was required to have px > 0.2 GeV and a 
charge opposite to that of the particle taken as a kaon. The cut was raised to 0.25 GeV 
for a data subsample, corresponding to an integrated luminosity of 16.9 ± 0.4 pb" 1 , for 
which the low-momentum track reconstruction efficiency was smaller due to the operating 
conditions of the CTD [30]. The corresponding D° candidate was assigned to a class of 
candidates "with AM tag" if the mass difference, AM = M(K7nr s ) —M(K7r), was in the 
range 0.143 < AM < 0.148 GeV. All remaining D° candidates were assigned to a class 
of candidates "without AM tag". For D° candidates with AM tag, the kaon and pion 
mass assignment was fixed by the track-charge requirements. For D° mesons without 
AM tag, the mass assignment is ambiguous. The pion and kaon masses can therefore be 
assigned to two tracks either correctly, producing a signal peak, or incorrectly, producing 
a wider reflected signal. To remove this reflection, the mass distribution, obtained for D° 
candidates with AM tag and an opposite mass assignment to the kaon and pion tracks, 
was subtracted from the M(Kn) distribution for all D° candidates without AM tag. 
The subtracted mass distribution was normalised to the ratio of numbers of D° mesons 
without and with AM tag obtained from a fit described below. 

Figure shows the M(K7r) distribution for D° candidates without AM tag, obtained 
after the reflection subtraction, and the M(Ktt) distribution for D° candidates with AM 
tag. Clear signals are seen at the nominal value of M(D°) in both distributions. The 
distributions were fitted simultaneously assuming the same shape for signals in both 
distributions. To describe the shape, a "modified" Gaussian function was used: 

Gauss mod oc exp[-0.5 • x^+W+o*-^ (1) 

where x = \ [M(Kir) — Mo]/a\. This functional form described both data and MC signals 
well. The signal position, M , and width, a, as well as the numbers of D° mesons in 
each signal were free parameters of the fit. Monte Carlo studies showed that background 
shapes in both distributions are compatible with being linear in the mass range above 
the signals. For smaller M(Ktt) values, the background shapes exhibit an exponential 
enhancement due to contributions from other D° decay modes and other D mesons. 
Therefore the background shape in the fit was described by the form [A + B • M(Ktt)] for 
M(Ktx) > 1.86 GeV and [A+B-M(KTr)]-exp{C-[M(K7r)-l.86]} for M(Kir) < 1.86 GeV. 
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The free parameters A, B and C were assumed to be independent for the two M(Kir) 
distributions. The numbers of D° mesons yielded by the fit were N untas (D°) = 11430±540 
and N tas (D°) = 3259 ± 91 for selections without and with AM tag, respectively. 

5.2 Reconstruction of additional D* + mesons 

The D* + D°7T+ events with p T (D* + ) > 3.8 GeV and \r](D* + )\ < 1.6 can be consid- 
ered as a sum of two subsamples: events with the D° having pt(D°) > 3.8 GeV and 
|?7(-D )| < 1.6, and events with the D° outside of that kinematic range. The former sam- 
ple is represented by D° mesons reconstructed with AM tag, as discussed in the previous 
section. The latter sample of "additional" D* + mesons was obtained using the same 
D° — > K~ix + decay channel and an independent selection described below. 

In each event, tracks with opposite charges and pt > 0.4 GeV were combined in pairs 
to form D° candidates. To calculate the pair invariant mass, M(Ktt), kaon and pion 
masses were assumed in turn for each track. Only D° candidates which satisfy 1.81 < 
M(Ktt) < 1.92 GeV were kept. Moreover, the D° candidates were required to have either 
p T (D°) < 3.8 GeV or \rj(D°)\ > 1.6. Any additional track, with p T > 0.2 GeV and a 
charge opposite to that of the kaon track, was assigned the pion mass and combined with 
the D° candidate to form a D* + candidate with invariant mass M(Ktttt s ). Here again the 
Pt cut was raised to 0.25 GeV for the data subsample for which the low-momentum track 
reconstruction efficiency was smaller. 

Figure 121 shows the AM distribution for the D* + candidates after all cuts. A clear signal 
is seen at the nominal value of M(D* + ) — M(D°). The combinatorial background was 
estimated from the mass-difference distribution for wrong-charge combinations, in which 
both tracks forming the D° candidate have the same charge and the third track has the 
opposite charge. The same tracks from a wrong-charge combination can produce two D° 
candidates due to an ambiguity in the kaon and pion mass assignment to tracks with the 
same charge. To exclude double counting, the multiple combinations of the same tracks 
which passed all cuts, including the M(Kir) requirement, were included with a weight 
1/2. 

The number of reconstructed additional D* + mesons was determined by subtracting the 
wrong-charge AM distribution after normalising it to the distribution of D* + candidates 
with the appropriate charges in the range 0.15 < AM < 0.17 GeV. The subtraction, 
performed in the signal range 0.143 < AM < 0.148 GeV, yielded A^ add (D* + ) = 826 ± 40. 

The AM distribution was also fitted to a sum of the modified Gaussian function (Eq. (^)) 
describing the signal and a threshold function describing the non-resonant background. 
The threshold function had a form A ■ (AM — m n ) B , where m n is the pion mass [8] 
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and A and B were free parameters. The results obtained using the fit instead of the 
subtraction procedure were used to estimate the systematic uncertainty of the signal 
extraction procedure. 

5.3 Reconstruction of D + mesons 

The D + mesons were reconstructed from the decay D + — > K~7t + tt + . In each event, 
two tracks with the same charges and pt > 0.5 GeV and a third track with opposite 
charge and px > 0.7 GeV were combined to form D + candidates. The pion masses were 
assigned to the two tracks with the same charges and the kaon mass was assigned to 
the third track, after which the candidate invariant mass, M^Kmr), was calculated. To 
suppress the combinatorial background, a cut of cos 9* (K) > —0.75 was imposed, where 
6*(K) is the angle between the kaon in the Ktttt rest frame and the Kirn line of flight 
in the laboratory frame. To suppress background from D* + decays, combinations with 
M(Kttit) — M(Kir) < 0.15 GeV were removed. The background from — > </>7r + with 
<f) — > K + K~ was suppressed by requiring that the invariant mass of any two D + candidate 
tracks with opposite charges was not within ±8 MeV of the <fi mass [4] when the kaon mass 
was assigned to both tracks. 

Figure El shows the M{Kttti) distribution for the D + candidates after all cuts. Reflections 
from and A+ decays to three charged particles were subtracted using the simulated 
reflection shapes normalised to the measured and A+ production rates. A clear signal 
is seen at the nominal value of D + mass. The mass distribution was fitted to a sum of a 
modified Gaussian function (Eq. (£[])) describing the signal and a linear function describing 
the non-resonant background. The number of reconstructed D + mesons yielded by the 
fit was N(D + ) = 8950 ± 600. 

5.4 Reconstruction of mesons 

The mesons were reconstructed from the decay Df — > 07r + with <p — > K + K~ . In each 
event, tracks with opposite charges and pr > 0.7 GeV were assigned the kaon mass and 
combined in pairs to form <p candidates. The <fi candidate was kept if its invariant mass, 
M(KK), was within ±8 MeV of the <fi mass [4]. Any additional track with px > 0.5 GeV 
was assigned the pion mass and combined with the <fi candidate to form a candidate 
with invariant mass M(KKtt). To suppress the combinatorial background, the following 
requirements were applied: 

• cosd* (it) < 0.85, where 0*(tt) is the angle between the pion in the KKtt rest frame 
and the KKtt line of flight in the laboratory frame; 
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• I cos 3 ^(K)! > 0.1, where 9'(K) is the angle between one of the kaons and the pion 
in the KK rest frame. The decay of the pseudoscalar meson to the <j) (vector) 
plus 7r + (pseudoscalar) final state results in an alignment of the spin of the <fi meson 
with respect to the direction of motion of the relative to D+. Consequently, the 
distribution of cos8'(K) follows a cos 2 6'(K) shape, implying a flat distribution for 
cos 3 ^(K). In contrast, the cosO'(K) distribution of the combinatorial background is 
flat and its cos 3 ^'^) distribution peaks at zero. The cut suppressed the background 
significantly while reducing the signal by 10%. 

Figure |U shows the M(KKn) distribution for the candidates after all cuts. Reflections 
from D + and A+ decays to three charged particles were subtracted using the simulated 
reflection shapes normalised to the measured D + and A+ production rates. A clear signal 
is seen at the nominal mass. There is also a smaller signal around the nominal D + 
mass as expected from the decay D + — > 07r + with — > K + K~. The mass distribution 
was fitted to a sum of two modified Gaussian functions (Eq. (JTJ)) describing the signals and 
an exponential function describing the non-resonant background. To reduce the number 
of free parameters, the width of the D + signal was constrained to 8/9 of the signal 
width; the constraint was verified by MC studies. The number of reconstructed Df 
mesons yielded by the fit was N(D+) = 1102 ± 83 5 . 

5.5 Reconstruction of A+ baryons 

The A+ baryons were reconstructed from the decay A+ — > K~pir + . In each event, two 
same-charge tracks and a third track with opposite charge were combined to form A+ 
candidates. Due to the large difference between the proton and pion masses, the proton 
momentum is typically larger than that of the pion. Therefore, the proton (pion) mass was 
assigned to those of the two tracks with the same charges which had larger (smaller) mo- 
mentum. The kaon mass was assigned to the third track and the candidate invariant mass, 
M(Kp-n), was calculated. Only candidates with pr(K) > 0.75 GeV, prip) > 1.3 GeV and 
Pt{^) > 0.5 GeV were kept. To suppress the combinatorial background, the following 
requirements, motivated by MC studies, were applied: 

• cos^*(A') > —0.9, where 9*{K) is the angle between the kaon in the Kptr rest frame 
and the Kpir line of flight in the laboratory frame; 

• cos#*(p) > —0.25, where 0*(p) is the angle between the proton in the Kpir rest frame 
and the Kpir line of flight in the laboratory frame; 

• p*(tt) > 90MeV, where p*(tt) is the pion momentum in the Kpir rest frame. 
5 The number of D + mesons, 239 ± 63, was not used further in the analysis. 
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To suppress the combinatorial background further, the measured dE/dx values of the 
three A+ candidate tracks were used. The parametrisations of the dE/dx expectation 
values and the Xi probabilities l p , Ik and of the proton, kaon and pion hypotheses, 
respectively, were obtained in the same way as described in a previous publication [23]. 
The l p , Ik and l n distributions for the A+ candidate tracks show sharp peaks around 
zero and become relatively flat towards one. To maximise the ratios of the numbers 
of correctly assigned protons, kaons and pions to the square roots of the numbers of 
background particles, the cuts l p > 0.15, Ik > 0.03 and l n > 0.01 were applied. The cuts 
rejected those ranges where the l p , Ik and l n distributions were at least twice as high as 
in the range 0.8 — 1. 

Figure 03 shows the M(Kpir) distribution for the A+ candidates after all cuts. Reflections 
from D + and decays to three charged particles were subtracted using the simulated 
reflection shapes normalised to the measured D + and production rates. A clear signal 
is seen at the nominal A+ mass. The mass distribution was fitted to a sum of a modified 
Gaussian function (Eq. describing the signal and a linear function describing the 
non-resonant background. The number of reconstructed A+ baryons yielded by the fit 
was A r (A+) = 1440 ± 220. 



6 Charm-hadron production cross sections 

The charm-hadron cross sections were calculated for the process ep — > e-D(A c )X in 
the kinematic region Q 2 < 1 GeV 2 , 130 < W < 300 GeV, p T (D,A c ) > 3.8 GeV and 
\r](D,A c )\ < 1.6. The cross section for a given charm hadron was calculated from 



A-C-B 



where N(D,A C ) is the number of reconstructed charm hadrons, A is the acceptance for 
this charm hadron, C is the integrated luminosity and B is the branching ratio or the 
product of the branching ratios [8] for the decay channel used in the reconstruction. The 
third uncertainties quoted below for the measured cross sections and charm fragmentation 
ratios and fractions are due to the branching-ratio uncertainties 6 . 

The combined Pythia+Rapgap MC sample was used to evaluate the nominal accep- 
tances. Small admixtures to the reconstructed signals from other decay modes were 
taken into account in the acceptance correction procedure. To correct from N ta,g (D°) 
(N untag (D )) to the production cross sections for D° mesons originating (not originating) 
from D* + decays, small migrations between the two samples were taken into account. 



6 Contributions from uncertainties of different branching ratios were added in quadrature. 
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The 6-quark relative contributions, predicted by the MC simulation using branching ra- 
tios of 6-quark decays to the charmed hadrons measured at LEP [31,32], were subtracted 
from all measured cross sections 7 . Subtraction of the 6-quark contribution reduced the 
measured cross sections by 3 — 7% and changed the measured charm fragmentation ratios 
and fractions by less than 4%. 

Using the reconstructed signals (see Section EJ) the following cross sections for the sum of 
each charm hadron and its antiparticle were calculated. The systematic uncertainties are 
discussed in Section |H1 

• the production cross section for D° mesons not originating from the D* + — > -D°7T+ 
decays: 

0- untag (Z>°) = 8.49 ± 0.44(stat.)t°;^(syst.)t°;20(br.) nb; 

• the production cross section for D° mesons originating from the D* + — > D 7i+ decays: 

a tas (D°) = 2.65 ± 0.08(stat.)t°;^(syst.) ± 0.06(br.) nb. 

The ratio a tas (D°)/Br)*+^D Tr+ gives the D* + cross section, a(D* + ), corresponding to 
D° production in the kinematic range pt{D°) > 3.8 GeV and \r](D°)\ < 1.6 for the 
D* + — > D°irf decay. Here Bd*+-*d°tt+ = 0.677 ±0.005 [8] is the branching ratio of the 
D* + -> D°7T+ decay; 

• the production cross section for additional D* + mesons: 

a add (D*+) = 1.05 ± 0.07(stat.)l^(syst.) ± 0.03(br.) nb. 

The sum a add (D* + ) + cr ta,g (D )/BD*+^£,o 7T + gives the production cross section for D* + 
mesons in the kinematic range pt{D* + ) > 3.8 GeV and |?7(-D* + )| < 1.6: 

a kil \D* + ) = 4.97 ± 0.14(stat.)l°;^( S y S t.)+0;}3(br.) nb; 

• the production cross section for D + mesons: 

a(D + ) = 5.07±0.36(stat.)+°;^(syst.)+°;^(br.)nb; 

• the production cross section for mesons: 

a(D+) = 2.37 ± 0.20(stat.) ± 0.20(syst.)+°;^(br.) nb; 

• the production cross section for A+ baryons: 

a(A+) = 3.59 ± 0.66(stat.)l°;^(syst.)iJ;^(br.) nb. 

7 The branching ratios of the 6-quark decays were updated using recent values [8] of the relevant charm- 
hadron decay branching ratios. 
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7 Charm fragmentation ratios and fractions 



7.1 Ratio of neutral to charged D-meson production rates 

Neglecting influences from decays of heavier excited D mesons, the ratio of neutral to 
charged .D-meson production rates is given by the ratio of the sum of D*° and direct D° 
production cross sections to the sum of D* + and direct D + production cross sections: 

a(D*°) + a dir (D°) 
u/d ~ a{D*+) + a^(D + Y 

where a dir (D°) and a dir (D + ) are those parts of the D° and D + inclusive cross sections 
which do not originate from D*° and D* + decays. Since all D*° decays produce a D° 
meson [8], the sum of a(D*°) and a dir (D°) is the production cross section for D° mesons 
not originating from D* + decays: 

a(D*°) + cr d[r (D°) = a unt ^(D°). (2) 

Subtracting from a(D + ) the contribution from D* + decays gives 

a dh (D+) = a{D+) - a(D* + ) ■ (1 - B D ^ D o w+ ). (3) 

Thus, the ratio of neutral and charged D-meson production rates can be calculated as 

p-untag ( Jj0^ ^untag ^ jjO\ 



Ru/d 



a(D+) + a(D*+) ■ B D *+^ D o T+ a(D+) + a^(D c 



Using the measured cross sections, the ratio of neutral to charged D-meson production 
rates, obtained for the kinematic region Q 2 < 1 GeV 2 , 130 < W < 300 GeV, pr(D) > 
3.8 GeV and \rj(D)\ < 1.6, is 

R u/d = 1.100 ± 0.078 (stat.)l°;°6! (syst.)iHlg (br.). 

The measured R u / d value agrees with unity, i.e. it is consistent with isospin invariance, 
which implies that u and d quarks are produced equally in charm fragmentation. 

Table Q compares the measurement with the values obtained in DIS [6] and in e + e~ 
annihilations. The latter value was calculated as 

R fjc -+ D°) - f(c -> D* + ) ■ B D «+^ D o n+ 

using fragmentation fractions compiled previously [7] and updated with the recent branch- 
ing ratio values [8] . All measurements agree with unity within experimental uncertainties. 
The branching ratio uncertainties of all measurements are highly correlated. 
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7.2 Equivalent phase-space treatment 



In the subtraction of the D* + contribution to D + production in Eq. (J3J), the cross-section 
a(D* + ), corresponding to D° production in the kinematic range pt{D°) > 3.8 GeV 
and \r](D°)\ < 1.6 for the D* + — > D°7rf decay, was used. Replacing a(D* + ) with 
a tag (D°) / B D *+^ D o n + gives 

a d "(D+) = a{D+) - a^(D°) ■ (1 - B D * + ^ + )/B D *+^ + . 

To compare direct D + and D* + production, the cross section a kin (D* + ) for pt(D* + ) > 
3.8 GeV and \r](D* + )\ < 1.6 is used in Section 17711 To compare the inclusive D + and D° 
cross sections with each other and with the inclusive D* + cross section it is necessary to 
take into account that only a fraction of the parent D* momentum is transfered to the 
daughter D meson. For such comparisons, the "equivalent" D + and D° cross sections 
were defined as the sums of their direct cross sections and contributions from D* decays 
calculated using (j kin (L)* + ) and ^(D* ): 

a c ^D + ) = a dil (D + ) + ( r kh \D*+)-(l-B D ^ D o w+ ), 
a C(l {D ) = a dil {D°) + a kin {D* + ) ■ B D *+^ D o n+ + a kin {D*°), 

where a km (D*°) is the inclusive D*° cross section for p T (D*°) > 3.8 GeV and |?7(_D* )| < 
1.6. This cross section can be written as the sum a(D*°) + a add (D*°), where a(D*°) is 
the part contributing to the D° production in the nominal kinematic range (as in Eq. (J2J)) 
and a add (D*°) is the production cross section for "additional" D*° mesons producing D° 
mesons outside of that kinematic range. The latter cross section was calculated using 
er add (-D* + ) and the expression for R u /d- 

untag/ r)0N 

a addp*0) = a addp*+) . Ru/d = a addp*+^ . 



a(D+) + a ta s(L) )' 

Using Eqs. (0) and for a dn (D°) and a dlT (D + ), respectively, and the expressions for 
a kiQ (D*°) and (T kin (L)* + ) gives 

a°*(D°) = a untag (D ) + ^(D°) + a add (D* + )-(i?^ + ^ + ^ + ), 
a cq ( J D+) = a( J D+) + ( 7 add ( J D*+).(l-^ + ^o 7r+ ). 

MC studies show that such "equivalent phase-space treatment" for the non-strange D and 
D* mesons minimises differences between the fragmentation ratios and fractions measured 
in the accepted Pt{D, A c ) and 77 (-D, A c ) kinematic region and those in the full phase space 
(see Section 7.6). 
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7.3 Strangeness-suppression factor 



The strangeness-suppression factor for charm mesons is given by the ratio of twice the 
production rate of charm-strange mesons to the production rate of non-strange charm 
mesons. All D* + and D*° decays produce either a D + or a D° meson, while all D* + 
decays produce a meson [8]. Thus, neglecting decays of heavier excited charm-strange 
mesons to non-strange charm mesons, the strangeness-suppression factor can be calculated 
as a ratio of twice the Df production cross section to the sum of D° and D + production 
cross sections. Using the equivalent D° and D + cross sections gives 

2a(Df) 2a(Df) 

Is 



a eq(D+) + a c i(D°) a(D+) + a xnskae (D°) + a ta s(D ) + <r add (D*+) • (1 + R u/d ) ' 

Using the measured cross sections, the strangeness-suppression factor, obtained for the 
kinematic region Q 2 < 1 GeV 2 , 130 < W < 300 GeV, p T {D) > 3.8 GeV and \rj(D)\ < 1.6, 
is 

7s = 0.257 ± 0.024 (stat.)lS (syst.)^ (br.). 

Thus, charm-strange meson production is suppressed by a factor w 3.9 in charm fragmen- 
tation. In simulations based on the Lund string fragmentation scheme [33], strangeness 
suppression is a free parameter which determines the ratio of probabilities to create s 
to u and d quarks during the fragmentation processes. In the absence of excited charm- 
strange meson decays to non-strange charm mesons, the Lund strangeness-suppression 
parameter would be effectively the observable, 7 S . In fact, production rates of the excited 
charm-strange mesons are poorly known; varying these rates in wide ranges in the PYTHIA 
simulation suggests that the Lund strangeness-suppression parameter is 10 — 30% larger 
than the observable, j s . 

Table 121 compares the measurement with the previous ZEUS 96-97 result, calculated from 
the ratio of Df to D* + cross sections [4], and with the values obtained for charm produc- 
tion in DIS [6] and in e + e~ annihilations. The e + e~ value was calculated as 

2f(c - Df) 

IS 



f(c^D+) + f(c^D°) 

using fragmentation fractions compiled previously [7] and updated with the recent branch- 
ing ratio values [8]. All measurements agree within experimental uncertainties. The large 
branching-ratio uncertainties are dominated by the common uncertainty of the — > 07r + 
branching ratio. This uncertainty can be ignored in the comparison with other measure- 
ments using the same branching ratios. 



13 



7.4 Fraction of charged D mesons produced in a vector state 

Neglecting influences from decays of heavier excited D mesons, the fraction of D mesons 
produced in a vector state is given by the ratio of vector to (vector+pseudoscalar) charm 
meson production cross sections. Only direct parts of the production cross sections for 
pseudoscalar charm mesons should be used. Using the expressions for a km (D* + ) and 
<7 dir (.D + ), the fraction for charged charm mesons is given by 

d = <t**(D*+) = a ta e(D )/B D . + ^ D o n+ +a add (D* + ) 

v a kinp*+) + o-dirp+j a(D+) + a l ^(D ) + a add (L>* + ) ' 

Using the measured cross sections, the fraction of charged D mesons produced in a vector 
state, obtained for the kinematic region Q 2 < 1 GeV 2 , 130 < W < 300 GeV, pr(D) > 
3.8 GeV and \rj(D)\ < 1.6, is 

P v d = 0.566 ± 0.025 (stat.)t°;°22 (syst.)^l (br.). 

The measured P d fraction is considerably smaller than the naive spin-counting prediction 
of 0.75. The predictions of the thermodynamical approach [34] and the string fragmenta- 
tion approach [35], which both predict 2/3 for the fraction, are closer to, but still above, 
the measured value. The BKL model [36,37], based on a tree-level perturbative QCD 
calculation with the subsequent hadronisation of the (c, q) state, predicts P d ~ 0.6 for 
charm production in e + e~ annihilations where only fragmentation diagrams contribute. 
For charm photoproduction, where both fragmentation and recombination diagrams con- 
tribute, the BKL prediction is P d ~ 0.66 in the measured kinematic range. 

Table El compares the measurement with the values obtained in DIS [6] and in e + e~ 
annihilations. The latter value was calculated as 

pd= f{c^D* + ) 

v /(c _ D+) + /(c ^ D * +) . B D * + ^ + 

using fragmentation fractions compiled previously [7] and updated with the recent branch- 
ing ratio values [8]. The measured P d value is smaller than, but consistent with, the pre- 
vious measurements. The branching-ratio uncertainties of all measurements are highly 
correlated. 

7.5 Charm fragmentation fractions 

The fraction of c quarks hadronising as a particular charm hadron, f(c—*D, A c ), is given 
by the ratio of the production cross section for the hadron to the sum of the production 
cross sections for all charm ground states that decay weakly. In addition to the measured 
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D°, D + , and A+ charm ground states, the production cross sections of the charm- 
strange baryons E+, 5° and should be included in the sum. The production rates 
for these baryons are expected to be much lower than that of the A+ due to strangeness 
suppression. The relative rates for the charm-strange baryons which decay weakly were 
estimated from the non-charm sector following the LEP procedure [38]. The measured 
S _ /A and tt~/A relative rates are (6.65 ± 0.28)% and (0.42 ± 0.07)%, respectively [8]. 
Assuming equal production of S° and S~ states and that a similar suppression is applicable 
to the charm baryons, the total rate for the three charm-strange baryons relative to the 
A+ state is expected to be about 14%. Therefore the A+ production cross section was 
scaled by the factor 1.14 in the sum of the production cross sections. An error of ±0.05 
was assigned to the scale factor when evaluating systematic uncertainties. 

Using the equivalent D° and D + cross sections, the sum of the production cross sections 
for all open-charm ground states (gs) is given by 

a gs = a^(D + ) + a c \D°) + a(D+) + *(A+) ■ 1.14, 

which can be expressed as 

a gs = a(D + ) + a untas (D°) + a tag (D°) + a add (/r+) • (1 + R u/d ) + a{Dj) + <r(A+) • 1.14. 

For the measured cross sections, 

a gs = 24.9 ± 1.0 (stat.)ii:I (syst.)t^ (br.) nb. 

The fragmentation fractions for the measured charm ground states are given by 

/(c-£>+) = a c \D + )/a gs = [a(D + )+a* dd (D* + )-(l-B D * + ^ D o n+ )]/a gs , 

f(c^D°) = a cq (D )/ a gs 

= [<t™***(D°) + <7 tag (-D°) + a* dd (D*+) ■ (R u/d + B D , + ^ + )]/a & , 

/(c-D+) = a{Dt)/a^ 

/(c-A+) = *(A+)/<v 

Using a km (D* + ), the fragmentation fraction for the D* + state is given by 

/(c /r + ) = a kin (D* + )/a gs = [a tag (D°) / B D *+^ D o 7T + + <7 add (_D* + )]/cr gs . 

The open-charm fragmentation fractions, measured in the kinematic region Q 2 < 1 GeV 2 , 
130 < W < 300 GeV, p T {D,A c ) > 3.8 GeV and |r/(D,A c )| < 1.6, are summarised in 
Table |U The results are compared with the values obtained in DIS [6] and with the 
combined fragmentation fractions for charm production in e + e~ annihilations compiled 
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previously [7] and updated with the recent branching-ratio values [8] . The branching-ratio 
uncertainties of all measurements are highly correlated. The measurements are consistent 
although the measured f(c — > D* + ) is smaller and /(c — > A+) is larger than those obtained 
in e + e~ annihilations. About half of the difference in the f(c — > D* + ) values is due to 
the difference in the /(c — » A+) values. The measurement may indicate an enhancement 
of A+ production in ep collisions with respect to e + e~ . However, this is unlikely to be 
a consequence of the baryon-number-flow effect [39] because no significant asymmetry 
between the A+ and A c production rates was observed 8 . 

7.6 Discussion of extrapolation effects 

The charm fragmentation ratios and fractions were measured in the region px{D,A c ) > 
3.8 GeV and \rj(D,A c )\ < 1.6. To minimise differences between the values measured in 
the accepted pr{D, A c ) and rj(D, A c ) kinematic region and those in the full phase space, 
the equivalent phase-space treatment for the non-strange D and D* mesons was used (see 
Section IT^j) . 

Table El shows estimates of extrapolation factors correcting the values measured in the 
accepted Pt(D, A c ) and f](D, A c ) region to the full phase space. The extrapolation factors 
were determined using three different fragmentation schemes: the Peterson parameteri- 
sation [40] of the charm fragmentation function as implemented in Pythia, the Bowler 
modification [20] of the LUND symmetric fragmentation function [21] as implemented in 
Pythia and the cluster model [22] as implemented in Herwig. The quoted uncertainties 
were obtained by varying relevant parameters in the Pythia and Herwig MC genera- 
tors. The extrapolation factors obtained are generally close to unity. The only exceptions 
are the factors given by the cluster model for f(c—> A+) and, to a lesser extent, for j s 
and /(c -£>+). 

This MC study suggests that the measured charm fragmentation ratios and fractions are 
close to those in the full Pt(D, A c ) and r\ (D, A c ) phase space. 

8 Systematic uncertainties 

The systematic uncertainties of the measured cross sections and fragmentation ratios and 
fractions were determined by changing the analysis procedure and repeating all calcula- 
tions. The following groups of the systematic uncertainty sources were considered: 

8 Separate fits of the M(K-pw + ) and M(K+pir-) distributions yielded N(A+)/N(A~) = 0.8 ± 0.2. 
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• {Si} the model dependence of the acceptance corrections was estimated using the 
Herwig MC sample, varying the Pt{D, A c ) and rj(D, A c ) distributions of the reference 
MC sample and by changing the MC fraction of charged D mesons produced in a vector 
state from 0.6 to 0.5 or 0.7; 

• {62} the uncertainty of the beauty subtraction was determined by varying the b- 
quark cross section by a factor of two in the reference MC sample and by varying the 
branching ratios of 6-quarks to charm hadrons by their uncertainties [31,32]; 

• {63} the uncertainty of the tracking simulation was obtained by varying all momenta 
by ±0.3% (magnetic field uncertainty), varying the track-loss probabilities by ±20% 
of their values and by changing the track momentum and angular resolutions by 

of their values. The asymmetric resolution variations were used since the MC signals 
typically had somewhat narrower widths than observed in the data; 

• {£4} the uncertainty of the CAL simulation was determined by varying the CAL energy 
scale by ±2%, by changing the CAL energy resolution by ±20% of its value and by 
varying the first-level trigger CAL efficiencies; 

• {6 5 } the uncertainties related to the signal extraction procedures were obtained as 
follows: 

— for the D° signals with and without AM tag: the background parametrisation 
and the range used for the signal fits were varied; 

— for the additional D* + signal: the range used for the background normalisation 
was varied or the fit was used instead of the subtraction procedure; 

— for the D + , and A+ signals: the background parametrisations, ranges used 
for the signal fits and amounts of the mutual reflections were varied. In addition, 
in the signal-extraction procedure, the constraint used for the D + — > KKn 
signal width was varied. In the A+ signal extraction procedure, an uncertainty in 
the dE/dx simulation was estimated by changing the dE/dx cut values in the MC 
and checking the effects with respect to changes expected from the \\ distribution. 

• {6 6 } the uncertainties of the luminosities of the e~p (±1.8%) and e + p (±2.25%) data 
samples were included taking into account their correlations; 

• {#7} the uncertainty in the rate of the charm-strange baryons (see Section I73J) . 

Contributions from the different systematic uncertainties were calculated and added in 
quadrature separately for positive and negative variations. The total and 61-67 systematic 
uncertainties for the charm-hadron cross sections and charm fragmentation ratios and 
fractions are summarised in Table El Correlated systematic uncertainties largely cancelled 
in the calculation of the fragmentation ratios and fractions. 



17 



To check the hadron-mass effects on the measured charm fragmentation ratios and frac- 
tions, the analysis was repeated using the charm-hadron transverse energy instead of the 
transverse momentum in the definition of the kinematic range of the measurement; the 
results obtained agreed with the reported values within statistical errors. The charm 
fragmentation ratios and fractions were also calculated separately for two W sub-ranges; 
no significant variations were observed. 

9 Summary 

The production of the charm hadrons D* + , D°, D + , and A+ has been measured 
with the ZEUS detector in the kinematic range pt(D,A c ) > 3.8 GeV, \r](D, A c )| < 1.6, 
130 < W < 300 GeV and Q 2 < 1 GeV 2 . The cross sections have been used to determine the 
charm fragmentation ratios and fractions with comparable precision to the e + e~ results. 

The ratio of neutral to charged D-meson production rates is 

R u/d = 1.100 ± 0.078 (statOlS (syst.)+°;^ (br.). 

The measured R u /d value agrees with unity, i.e. it is consistent with isospin invariance, 
which implies that u and d quarks are produced equally in charm fragmentation. 

The strangeness-suppression factor is 

7s = 0.257 ± 0.024 (stat.)lffi (syst.)iffi. ( br 0- 

Thus, Ds-meson production is suppressed by a factor ps 3.9 in charm fragmentation. 
The fraction of charged D mesons produced in a vector state is 

P v d = 0.566 ± 0.025 (stat.)iffi (syst.)!™ 2 , 2 (br.). 

The measured fraction is considerably smaller than the naive spin-counting prediction of 
0.75. The predictions of the thermodynamical approach [34] and the string fragmenta- 
tion approach [35], which both predict 2/3 for the fraction, and the BKL model [36,37] 
prediction (« 0.66) are closer to, but still above, the measured value. 

The measured R u /d and 7 S values agree with those obtained in DIS [6] and in e + e~ 
annihilations. The e + e~ values were calculated using fragmentation fractions compiled 
previously [7] and updated with the recent branching ratio values [8]. The measured 
value is smaller than, but consistent with, the previous measurements. 

The fractions of c quarks hadronising as D* + , D°, D + , and hadrons have been 
calculated in the accepted kinematic range. The measured open-charm fragmentation 
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fractions are consistent with previous results although the measured f(c — > D* + ) is smaller 
and f(c — > A+) is larger than those obtained in e + e~ annihilations. About half of the 
difference in the f(c — > D* + ) values is due to the difference in the f(c — > A+) values. 

These measurements generally support the hypothesis that fragmentation proceeds inde- 
pendently of the hard sub-process. 
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R u /d 


ZEUS (tp) 


1.100 ± 0.078(stat.)^;^(syst.)^9(br.) 


HI (DIS) [6] 


1.26 ± 0.20(stat.) ± 0.11(syst.) ± 0.04(br. © theory) 


combined e + e~ data [7] 


1.020 ± 0.069(stat. © syst.)+°;°^(br.) 



Table 1: The ratio of neutral to charged D-meson production rates, Ru/d- 





7s 


ZEUS ( 7 p) 


0.257 ± 0.024(stat.)i°;^(syst.)lHl9(br.) 


ZEUS 96-97 [4] 


0.27 ± 0.04(stat.)+°;^( sys t.) ± 0.07(br.) 


HI (DIS) [6] 


0.36 ± 0.10(stat.) ± 0.01(syst.) ± 0.08(br. © theory) 


combined e + e~ data [7] 


0.259 ± 0.023(stat. © syst.)+°;°^(br.) 



Table 2: The strangeness- suppression factor in charm fragmentation, 7 S . 





pd 

V 


ZEUS ( 7 p) 


0.566 ± 0.025(stat.)™(syst.)lo:o23(br.) 


HI (DIS) [6] 


0.693 ± 0.045(stat.) ± 0.004(syst.) ± 0.009(br. © theory) 


combined e + e~ data [7] 


0.614 ± 0.019(stat. © syst.)+°;°|(br.) 



Table 3: The fraction of charged D mesons produced in a vector state, P^. 
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ZEUS hp) 
p T (D,A c ) > 3.8 GeV 
\V(D,A C )\ < 1.6 


Combined 
e + e~ data [7] 


HI (DIS) 
[6] 




stat. syst. br. 


stat.© syst. br. 


total 


/(c ->£>+) 


n 217 ± n ni4 +o-° l3 +o-o 14 

U.Z1I -0.005-0.016 


0.226 ±0.010 t°o° u 


0.203 ±0.026 


/(c-D°) 


0.523 ±0.021 i^l™ 2 . 2 . 


0.557 ±0.023 


0.560 ±0.046 




0.095 ±0.008 ^roTr 


0.101 ±0.009 tffi 


0.151 ±0.055 


/(c-A+) 


0.144 ±0.022 ±g:SS±85S 


0.076 ±0.007 tffi 




f(c^D*+) 


0.200 ±0.009 t° Zt° Z 


0.238 ±0.007 ±g:gg| 


0.263 ±0.032 



Table 4: The fractions of c quarks hadronising as a particular charm hadron, 
/(c — * D, A c ). TTie fractions are shown for the D + , D° , and A + charm ground 
states and for the D* + state. 





Peterson 


Bowler 


Cluster model 




(Pythia) 


(Pythia) 


(Herwig) 


Ru/d 


n qq+0.02 
u - yj -o.oo 


n qq+ - 02 

u - JJ -0.00 


i,uu -0.00 


Is 


1 04 +0 - 04 


1 oo +0 - 05 

i - uu -0.04 


i i O+0.07 
i - L °-0.05 


pd 

V 


1.00 ±0.02 


q7 +0 - 01 
U - J ' -0.00 


q6 +a02 
u - JU -o.oi 


f(c^D+) 


1 oo +0 - 02 


i 02++ - 01 


n qq+ 01 


f(c - D°) 


0.99 ±0.01 


0.98 ±0.01 


qfi+ 00 

u - yu -0.02 




1 m+ - 03 

± - uo -0.06 


1 oo +0 - 04 

i - uu -0.03 


i i c+0.06 
i - id -0.05 


/(c-A+) 


1 m+ - 02 

i - ui -0.05 


i n8+ - 03 




/(c^L>*+) 


1 00 +002 

i - uu -0.03 


qfi+ - 00 

u - JU -0.02 


q3+ 01 

u - yo -0.02 



Table 5: The estimates of extrapolation factors which correct charm fragmentation 
ratios and fractions measured in the accepted Pt{D , A c ) andrj(D,A c ) region to the 
full phase space. For further details, see text. 
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total 


5l 

w 1 


$2 


5 3 

^ o 


5a 


<$5 

w 


s e 

v o 


5 7 




(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


untag/ r)0\ 


+5.5 
-5.6 


+2.8 
-0.6 


+1.8 
-3.4 


+1.1 

-1.4 


+1.3 
-1.2 


+3.4 
-3.4 


+2.2 
-2.2 




a yu ) 


+4.0 
-3.7 


+2.5 
-1.2 


+1.1 
-2.1 


+1.4 
-1.3 


+1.4 
-1.1 


+0.7 
-0.4 


+2.2 
-2.2 




u yi^> ) 


+8.4 
-3.6 


+5.8 
-0.4 


+1.0 
-1.9 


+3.3 
-1.6 


+1.8 
-1.4 


+4.2 
-0.1 


+2.2 
-2.2 






+4.6 
-3.6 


+3.1 
-1.0 


+1.1 
-2.0 


+1.6 
-1.3 


+1.4 
-1.1 


+1.1 
-0.3 


+2.2 
-2.2 




a(D±) 


+8.7 
-4.5 


+4.5 
-0.3 


+1.8 
-3.3 


+1.0 
-1.6 


+1.3 
-1.0 


+6.7 
-0.6 


+2.3 
-2.3 






+8.3 
-8.5 


+6.0 
-0.0 


+3.9 
-7.0 


+1.4 
-1.3 


+1.8 
-1.0 


+2.9 
-4.0 


+2.2 
-2.2 






+15.1 
-18.3 


+13.2 
-0.6 


+3.1 
-5.6 


+4.8 
-1.5 


+2.2 
-3.0 


+3.3 
-16.9 


+2.3 
-2.3 




R u /d 


+3.5 
-5.5 


+0.0 
-0.9 


+0.4 
-0.7 


+0.6 
-0.6 


+0.1 
-0.3 


+3.4 
-5.4 


+0.2 
-0.2 




Is 


+5.0 
-6.3 


+2.4 
-0.2 


+2.4 
-4.1 


+0.9 
-0.8 


+0.8 
-0.0 


+3.5 
-4.7 


+0.1 
-0.1 




pd 

V 


+1.2 
-3.9 


+0.4 
-1.0 


+0.8 
-0.5 


+0.5 
-0.1 


+0.1 
-0.1 


+0.6 
-3.7 


+0.1 
-0.1 






+6.8 
-5.7 


+5.3 
-0.4 


+2.0 
-3.8 


+1.8 
-1.0 


+1.2 
-1.1 


+1.9 
-3.1 


+2.2 
-2.2 


+0.7 
-0.7 


f(c^D+) 


+6.1 
-2.1 


+0.3 
-0.8 


+0.7 
-0.4 


+0.2 
-1.0 


+0.6 
-0.3 


+6.0 
-1.5 


+0.1 
-0.1 


+0.7 
-0.7 


f(c -> D°) 


+3.4 
-3.2 


+0.2 
-2.2 


+0.9 
-0.5 


+0.3 
-0.8 


+0.5 
-0.4 


+3.1 
-2.1 


+0.1 
-0.1 


+0.7 
-0.7 




+4.9 
-5.4 


+0.8 
-0.2 


+1.9 
-3.3 


+0.3 
-0.8 


+1.2 
-0.3 


+4.1 
-4.1 


+0.1 
-0.1 


+0.7 
-0.7 


/(c - A+) 


+9.1 
-15.1 


+7.4 
-0.3 


+1.2 
-1.9 


+3.5 
-0.7 


+1.7 
-2.9 


+3.1 
-14.6 


+0.2 
-0.2 


+0.7 
-0.7 


f(c -> 


+3.9 
-3.2 


+0.4 
-2.2 


+1.9 
-1.0 


+0.6 
-0.5 


+0.7 
-0.5 


+3.1 
-1.9 


+0.1 
-0.1 


+0.7 
-0.7 



Table 6: The total and Si-8j (see text) systematic uncertainties for the charm 
hadron cross sections and charm fragmentation ratios and fractions. 
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Figure 1: The M(Kir) distributions (dots) for (a) the D°/D° candidates without 
AM tag, obtained after the reflection subtraction (see text), and for (b) the D°/D° 
candidates with AM tag. The solid curves represent a fit to the sum of a modified 
Gaussian function and a background function. 
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Figure 2: The distribution of the mass difference, AM = M(Kttit s ) — M(Ktt) , for 
the "additional" D*^ candidates (dots). The histogram shows the AM distribution 
for wrong-charge combinations. The shaded band shows the signal range in which 
the wrong-charge background subtraction was performed. 
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Figure 3: The M{Kttti) distribution for the D candidates (dots). The solid 
curve represents a fit to the sum of a modified Gaussian function and a linear 
background function. 
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Figure 4: The M(KKtt) distribution for the Df candidates (dots). The solid 
curve represents a fit to the sum of two modified Gaussian functions and an expo- 
nential background function. 
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Figure 5: The M(Kpir) distribution for the A+/A c candidates (dots). The 
solid curve represents a fit to the sum of a modified Gaussian function and a linear 
background function. 
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